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Abstract 

The magnetic molecule Fes has been observed to have a rich pattern of 
degeneracies in its magnetic spectrum as the static magnetic field applied to 
the molecule is varied. The points of degeneracy, or diabolical points in the 
magnetic field space, are found exactly in the simplest model Hamiltonian for 
this molecule. The points are shown to form a perfect centered rectangular 
lattice, and are shown to be multiply diabolical in general. The multiplicity 
is found. An earlier semiclassical solution to this problem is thereby shown 
to be exact in leading order in 1/J where J is the spin. 
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A famous theorem of quantum mechanics states that the intersection of two energy levels 
of a physical system is infinitely unlikely as a single parameter is varied [|l|]. Instead, level 
repulsion is the rule, and level crossing happens only when there is some symmetry. If 
one can vary more than one parameter, however (at least two if the Hamiltonian is real, 
three if it is complex), then degeneracies can be found at isolated points in the parameter 
space Such degeneracies are usually called "accidental", for as a rule the points obey 

no discernible pattern. Berry and Wilkinson's coinage for such degeneracies, "diabolical," 
emphasizes the same point, although the term originates less in any desire to attribute the 
phenomenon to a Mephistophelean hand, than in the fact that the energy surface near the 
degeneracy is a double cone, like an Italian toy called the diavolo. When exceptions occur, 
and the pattern of diabolical points or degeneracies is not random, as in the Kepler problem, 
or in the isotropic harmonic oscillator, one seeks (and sometimes finds) a higher dynamical 
symmetry of the Hamiltonian. 

It is therefore of some interest that a whole array of diabolical points has been discov- 
ered by Wernsdorfer and Sessoli (WS) in the magnetic spectrum of the molecular solid 
[(tacn^FesC^OH)^] 8-1- (shortened to Fes from now on). Strong intramolecular exchange 
interactions among the Fe 3+ ions lead to a total spin of J = 10 in the ground manifold for 
each molecule, which can therefore be conceived of as a single large spin of this magnitude. 
The dipolar interactions between molecules are weak, and may be ignored. Spin-orbit inter- 
actions lead to an easy-axis anisotropy, with further differentiation in the non-easy plane. 
The simplest anisotropy Hamiltonian that describes the magnetic properties of one molecule 
is 



with h\ > k2 > ( ki pa 0.33 K, and k% ~ 0.22 K §), g = 2, and H is an external magnetic 
field. With Eq. ([I]), the easy, medium, and hard axes are z, y, and x, respectively. 

To understand the evidence for diabolicity, let us suppose that the field H is along the 
hard axis, x, and not too large. The easy directions, which are along ±z when H = 0, 
cant symmetrically toward x. At first sight, as H x is increased, it should be progressively 
easier for the spin to tunnel from a state localized in the potential well in the positive J z 
hemisphere to the symmetrically located state in the negative hemisphere. It was found 
some time ago, however, that in the model (|l|), the tunnel splitting between the ground 
states in the positive and negative J z wells actually oscillates as a function of H x , going 
to zero at periodically spaced H x values 0. Precisely these oscillations have been seen by 
WS in Fes- The oscillations were initially explained in terms of an interference between 
instanton trajectories for the spin, and the observations on Fes are very hard to explain 
by other mechanisms. The clincher is that WS see additional oscillations that were not 
predicted. To understand these, suppose H also has a nonzero z component so as to bring 
the first or second excited state in the positive J z well into degeneracy with the ground state 
in the other well. One can then conceive of tunneling between these states. WS observe 
that the amplitude for this tunneling also oscillates with H x , and that the oscillations are 
shifted by half a period for each excited state in the deeper well. 

That the tunnel splitting between two states vanishes, is merely another way of saying 
that the states are exactly degenerate. The degeneracies that lie on the H x or H z axes in 
the magnetic field space can be understood in terms of symmetry allowed level crossings as 




(1) 
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per the von Neumann- Wigner theorem , but the newly discovered ones by WS lie off the 



axes, and cannot be so understood. They are truly nontrivial instances of diabolical points. 
Since the experiments were reported, they have been successfully explained by one of us 
(AG) [0, and independently, Villain and Fort fll2| . The approach is based on a discrete 



version of the phase integral or Wentzel-Kramers-Brillouin method, and is semiclassical in 
nature, with 1/J playing the role of h in continuum problems with massive particles. It is 
found that the £'th level in the negative J z well (with £' = being the lowest level) and the 
th level in the positive one are degenerate when H y = 0, and (see Fig. |l|) 

H z (£',£") V\{£"-£')_ (2) 
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with n = 0,1, ... ,2J -(£' + £" + 1). Here, A = k 2 /h, and H c = 2k x Jjg\x B . Note that these 
equations do give a half-period shift per excited state, as seen by WS. 

The surprise is that although Eqs. (Q) and ([3]) are dervied semiclassically, and should 
have higher order corrections in 1/J, they appear to be exact as written! This has been 
noted by both Villain and Fort, and AG. The evidence is from (a) analytic diagonalization 
for J < 2, (b) perturbation theory in A [13], and (c) numerics. Note that if exact, Eqs. (0) 



and would imply not only that the diabolical points lie on a perfect centered rectangular 
lattice in the H x -H z plane, but also that many of the points are multiply diabolical, i.e., 
that more than one pair of levels is simultaneously degenerate. It is easily shown that the 
multiplicity is as indicated in Fig. [IJ If we arrange the points into concentric rhombi, those 
on the outermost rhombus are singly diabolical (i.e., there is only one pair of degenerate 
states), those on the next rhombus are doubly diabolical (two pairs of degenerate states), 
and so on. 

In this paper we shall prove that this perfect lattice hypothesis is in fact true. Not only 
is this an interesting problem in mathematical physics in its own right, but we believe that it 
will help understand real Fe 8 also. Exact solutions generally open the way for perturbative 
treatment of small corrections, and in Fes, we beleive our work will enable us to better treat 
such effects as the higher order anisotropies and the dipolar interactions mentioned above, 
and allow more detailed understanding of nonzero temperature effects [[TJj] . 

Before giving our proof, we should note that we have not been able to find if the Hamil- 
tonian ([IJ) has a higher symmetry at the diabolical points. Such a suspicion is natural, 
given the experience with the exceptional cases mentioned in our first paragraph. Likewise, 
we have not succeeded in finding the wavef unctions. We also note that the semiclassical 
approximation is demonstrably inexact at non diabolical values of H, in order to dispel any 
suspicion in the readers' minds that the classical and quantum dynamics of the Hamiltonian 
(|I|) are identical as in the case of the harmonic oscillator. 

We now present our proof. It proceeds in three steps. In step 1 we perform a spin rotation 
about y so that 7i no longer has any terms in Jj. The Hamiltonian is then tridiagonal 
in the new J z basis. In step 2 we make use of a necessary condition for a tridiagonal 
Hermitean matrix to have degenerate eigenvalues, and thus determine the possible locations 
of any diabolical points. In step 3, we use continuity and topological arguments to find the 
multiplicity of each of these diabolical points. The results are precisely those given above. 
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Step 1: Let us rotate J about the y axis so that 




Vx Vx\(j z \ 
-Vx Vx ) \Jx) ' 



(4) 



where A = 1 — A. It is also convenient to define scaled fields u x , u y , and u z via 



H = ^j(V\u x , u y , s/Xu x 



(5) 



to scale all energies by ki, and to write TC = TL/kx- In the new axes, we have 
H = -(\- \)J 2 Z - \f)^(J z J x + J X J Z ) 

- (XU z -\U x )J z + \[}J\{U X + U Z )J X +UyJy 



(6) 



Step 2: The Hamiltonian (^) is tridiagonal in the the new J z basis, J z \m) = m\m). Its 
only nonzero matrix elements are (m\H\m) = w m , and (m\H\m!) = t m ^ m i, with mf = mil. 
And, it is obviously Hermitean: = w m , t m>m i = t* m , m . For such matrices, it is a theorem 
that all the eigenvalues are simple (i.e., nondegenerate) if none of the t m ^ m i vanish |L5 |. This 



result is physically almost obvious if one thinks of Tt as a tight-binding model for an electron 
in one dimension with on-site energies w m and nearest neighbor hopping elements £ m ,m±i- 
Two states which were degenerate could be spatially localized in different regions of the 
lattice. But then, since t m , m ±i ^ for any m, it would be possible for the electron to hop 
from one region to the other, which is self-contradictory. The rigorous proof by Stoer and 
Bulirsch consists of noting that the successive diagonal subdeterminants of H — EI (where 
I is the unit matrix), form a Sturmian sequence of polynomials pj(E), j = 1, 2, . . . , 2 J + 1, 
with the properties that pj is of degree j, has j real roots, each one of which is simple, and 
is strictly bracketed by two roots of Pj+i- 

It follows that the diabolical points of 7i, if any, must lie on the loci in magnetic field 
space defined by t m)m+ i = 0. Using the standard representation of the angular momentum 



matrices, we have 



t 



m,m+l 



[u x + u z + (2m + 1)) v AA + iu y 



[J{J 



m(m + 1)] 1/2 . 



(7) 



The real and imaginary parts of this quantity must vanish separately for some m. We thus 
conclude that any diabolical points must lie in the H x -H z plane: 



u, 



0. 



and in this plane on the lines 



u x + u z = — (2m + 1), m = — J, — J + 1, . . . , J — 1. 



(8) 



(9) 



Further, when these conditions are obeyed, TC divides into two blocks, of size n n and n p , 
with n n = J + m + 1, and n p = J — m . In the first block, m < m , and in the second 
m > m Q + 1. 
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The conditions @ and (|9|) are not enough to locate the diabolical points. However, we 
can repeat our argument, starting with a rotation (|J) in which the signs of the Vx entries 
are reversed. In this way we find that at a diabolical point, we must also obey the condition 

u x - u z = -(2n + 1), n = -J,-J + l,...,J-l, (10) 

which along with Eq. @ fixes the points completely. Solving these equations, we get 

u x = -{m + n + 1), 
u z = n - wo, 

which with Eq. (||) are identical to Eqs. (H) and (|3|). 

Step 3: It remains to find just how many pairs of levels are degenerate at each of the 
points (fLip . To do this, we first show that under a continuous change in the parameter A, 
the diabolical points of the Hamiltonian (|l|) must evolve smoothly in the u x -u z plane. In 
particular, the number of diabolical points must stay fixed. We present two arguments for 
this, one more physical, and the other more mathematical. 

The physical argument is much like the standard one for the von Neumann- Wigner 
theorem. Suppose that for some A = Aq, Ti has two degenerate states \if) a ) and \ipb) at 
Ux = u x o, u z = u z q. (We set u y = throughout.) Now let A be changed by a small amount 
SX. We seek new eigenstates in the first approximation as linear combinations of \ip a ) and 
\ipb)- The secular matrix in this subspace is given by 



7~tba 7~Lbb 



(12) 





with Hba = 7~tab since with u y = 0, the Hamiltonian is real. Here, H aa , 7~Lab, etc., are all 
smooth functions of A, u x , and u z . The eigenvalues of this matrix will be equal if and only 
if 'Haa = 7~Lbb, and H a b = 0. Expanding these conditions in the deviations from A , u x0 , and 
u z0 , we get 

(13) 

where the quantities a through d are partial derivatives of TL aa etc. with respect to u x and u z , 
and e and / are partial derivatives with respect to A. Since the diabolical point for A = Ao is 
isolated, it follows that the only solution to Eq. (|TjD when 5X = is 5u x = 5u z = 0, i.e., that 
the matrix on the left hand side is nonsingular, and a nonzero solution can be found when 
5X 7^ 0. It follows in turn that the diabolical point can only be moved, and not eliminated 
or created by changing A. 

The second argument is based on Berry's phase |Tj|. Let us consider a diabolical point as 
in the preceding paragraph, and let C be a small closed contour in the u x -u z plane encircling 
this point. Berry's phase is given by 

7 (C) = i I (Va(A,u)|V u ^ a (A,u)> ■ du, (14) 

where u = (u x ,u z ), and V u is a gradient with respect to these fields. As shown by Berry, 
7(C) = ±7r if C encloses a true diabolical point, and 7(C) = if the two states merely 
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approach each other very closely without ever being degenerate. [Actually, since our Hamil- 
tonian is real, and the parameter space (u x ,u z ) is two-dimensional, we really only need the 
weaker result due to Herzberg and Longuet-Higgins H for the sign change of the wavefunc- 
tion upon encircling the degeneracy: e* 7 ^ = —1.] Since the perturbation 87i engendered 
by changing A or u by a small amount is non-singular, \4> a {\ u)) is a smooth function of 
A and u. It follows that if A varies continuously, the integrand of Eq. ([14]) can not change 
discontinuously. Hence, for small enough 5X, the phase 7(C) must continue to be what it 
was for A = A , +tt, say, implying that C continues to encircle a degeneracy. 

The rest is plain sailing. Consider 7i in the form (0), and let u x and u z be one of the 



points (|TT|). As noted before, 7i divides into two blocks at this point. Consider now how 
the eigenvalues of each block change as A is varied. Suppose the variation is as depicted in 
Fig. 0. This would imply that the total diabolicity of our system changes discontinuosly, 
in violation of the result just proved. Thus a behavior as in Fig. ^| is impossible, and the 
correct picture is as in Fig. |3|. One can also see that the behavior can not be that as in 
Fig. ^| for some (mo, no) and that as in Fig. ^| for others, since there is then no way to keep 
the total diabolicity constant. 

We thus conclude that the multiplicity f(m ,n ) of the diabolical point at (m ,n ) is 
independent of A, and we may find it by evaluating it for any one value of A. We do this for 
A = 0, as Ti. is then not just tridiagonal, but diagonal. Degeneracy occurs whenever 

w mi =w m2 , (mi < m ; m 2 > m + 1). (15) 

Since w m = m 2 — (m — n Q )m for A = 0, this condition becomes 

mi + m 2 = m — no- (16) 

The problem is now one of counting. We leave it as an exercise to show that with the 
conditions in Eq. (|l~5| ) on mi and m-2, and in Eqs. (|D|) and ([K]) on itlq and no, the number 
of solutions to Eq. (|i~6|) is given by 



f(m ,n ) = -[2 J + 1 - |m - n Q \ - \m + n + 1|]. (17) 

This is precisely the multiplicity implied by Eqs. (||) and (||). One way to see this is to 
evaluate / on the rhombi drawn in Fig. p], in particular, the segment lying in the first 
quadrant in the H x -H z plane. Regarding the outermost rhombus as the first, the next as 
the second, etc., we have m = — ( J + 1) + k on the kth one. Further, for the part in the 
first quadrant, n ranges from m to J — k. Thus (m — n ) and (m + n + 1) are both 
negative, and 

/(m ,n Q ) = ^[2J + 1 + (m - n ) + (m + n + 1)] 
= J + m + 1 

= k, (18) 

exactly as expected. 

For completenes, we conclude by observing that our results obey the symmetries of the 
Hamiltonian. First, for half-integral J, the point at H = corresponds to mo = no = —1/2. 
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The value of / is then J+ ~, i.e., every state is doubly degenerate, as required by Kramers' 
theorem. Second, the reflection symmetries H x — > —H x , H z — > —H z are clearly obeyed by 
the set (0) and (|3|). In terms of the quantities mo and uq, these symmetries correspond to 



mo — ► —(mo + 1) and n — > —(no + 1), and it is easy to see from Eq. (|17D that these leave 
/ unchanged. Lastly, the Hamiltonian ([[]) has the following duality property. Showing its 
dependence on A, H x and H z explicitly by writing 7i(A, H x , H z ), we note that a 90° rotation 
about (x + z)/ yields the transformation 

H(\, H x , H z ) ^-H(l- A, H„ H x ). (19) 

In particular, the spectra of the two Hamiltonians are so related, and ranking the levels is 
order of increasing energy, we see that if the levels with ordinal numbers i and i + 1 are 
degenerate when H x = f x (X) and H y = f y (X), where the functions f x and f y are unknown, 
then level numbers 2 J + 2 — i and 2 J + 1 — i are degenerate when H x = f y (l — A) and 
H y = f x {l — X). This property is also obeyed by Eqs. (||) and @, and (|l7j). 
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FIGURES 



FIG. 1. Diabolical points for the model Hamiltonian (|l|), shown for the case J = 7/2. The 
points can be organized into rhombi (including the degenerate central rhombus with only one point). 
All points on a rhombus have the same multiplicity, i.e., number of simultaneously degenerate levels. 
These numbers are as shown. The dashed lines are self-dual under the transformation (|l~9|). 

FIG. 2. Impermissible variation of the energy levels with A for fixed rriQ and no- The solid 
curves are the energies of one of the blocks of TL, and the dashed curves of the other. 

FIG. 3. Same as Fig. ||, except that the variation is now permissible. The solid and dashed 
curves drawn right next to each other are in fact coincident. 
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